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Abstract 

The thorough analysis of the ALEPH data Q on hadronic r-decay is performed in 
the framework of QCD. The perturbative calculations are performed in 3 and 4-loop 
approximations. The terms of the operator product expansion (OPE) are accounted up 
to dimension D = 8. The value of the QCD coupling constant as(m^) = 0.355 it 0.025 
was found from hadronic branching ratio Rr- The V + A and V spectral function are 
analyzed using analytical properties of polarization operators in the whole complex g^- 
plane. Borel sum rules in the complex q'^ plane along the rays, starting from the origin, 
are used. It was demonstrated that QCD with OPE terms is in agreement with the data 
for the coupling constant close to the lower error edge a<j(m^) = 0.330. The restriction 
on the value of the gluonic condensate was found i^G^) = 0.006 ± 0.012 GeV^. The 
analytical perturbative QCD was compared with the data. It is demonstrated to be in 
strong contradiction with experiment. The restrictions on the renormalon contribution 
were found. The instanton contributions to the polarization operator are analyzed in 
various sum rules. In Borel transformation they appear to be small, but not in spectral 
moments sum rules. 



PACS: 13.35.D, 11.55.H, 12.38 



1 Introduction 



The high precision data on hadronic r-decay, obtained by ALEPH [|l] , OPAL and CLEO 
collaborations, namely the measurements of the total hadronic branching ratio Rr = 
B{t Ut- + hadrons)/-B(r eVe^r)) vector V and axial A spectral functions allow one to 
perform various tests of QCD at low energies: to determine asiQ"^) at low Q^, to check the 
operator product expansion (OPE) and to perform search for other possible nonperturbative 
modifications of QCD — renormalons, analytical as(Q^), instantons etc. An early attempt to 
check OPE in QCD based on e^e~ annihilation data has been made by Eidelman, Vainstein 
and Kurdadze 0] but the accuracy of the data at that time was not good enough. Also the 
authors of took as granted that the QCD coupling constant is rather small, Ag^^^ (for 3 
flavors) is about 100 MeV and neglected higher order terms of perturbative series. Now it is 
common belief that is much larger and A'-^-' ~ 300 — 400 MeV in 2-3 loop approximation. 
Therefore the problem deserves reconsideration. 

In the previous paper by two of us (B.I. and K.Z.) the difference of vector and axial 
current correlators was analyzed using ALEPH data on r-decay |Q . The analytical properties 
of the polarization operator in the whole complex g^-plane were exploited and the vacuum 
expectation values of dimension 6 and 8 operators (vacuum condensates) were found. Here 
we consider V + A correlator, where perturbative corrections are dominant. 

Deflne the polarization operators of hadronic currents: 

KM = ^ I e^"" (TJ^(x) J.(O)t) dx = {q,q. - g,.q') H«(g2) + q,q.Uf{q') , (1) 

where J = V,A; V^ = u-f^d , A^ = u-ff^-f^d . 
The imaginary parts of the correlators are the so-called spectral functions [s = g^), 

vi/ai{s) = 27TlmU'yj^{s + iO) , ao{s) = 2n ImUf {s + iO) . (2) 

which have been measured from hadronic r-decays for < s < m^. 

The spin-1 parts Ily\q'^) and H^^(g^) are analytical functions in the complex g^-plane 
with a cut along the right semiaxes starting from the threshold of the lowest hadronic state: 
4m^ for Hy"* and 9m^ for Ka ■ The latter has a kinematical pole at = 0. This is a specific 
feature of QCD, which follows from the chiral symmetry in the limit of massless u, rf-quarks 
and its spontaneous violation. It can be easily shown (see also [^), that the kinematical 
pole arises from the pion contribution to H^^, which is given by 




where /^r is the pion decay constant, f-^ = 130.7 MeV 



2 Hadronic branching ratio and the value of as{rri^) 

The total hadronic branching ratio into final state with zero strangeness is given by well 
known expression, which can be written in the following form (see e.g. [§]): 

B{t ^ Ut- + hadrons5=o) 
Rt,v+a = ^7 
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where \Vud\ = 0.9735 ± 0.0008 ||^ is Cabbibo-Kabayashi-Maskava matrix element, Sew = 
1.0194±0.0040 includes electroweak corrections 0. The spin-0 axial spectral function a^ls) is 
basically saturated by r — >• ttz/^ channel and can be read off from (^): ao{s) = 27r^/^5(s— m^). 
So the last term in (HI) gives small correction 
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n^''(s), for which the short 



The rest of contains only the imaginary part of 11 

notation n(s) will be used later on. As follows from U^^^q'^) compensates the kinematical 
pole at = in n^'*(g^). So the combination n(g^) has no kinematical poles and is an 
analytical function of in the complex g^-plane with a cut along the positive real axis. 

The convenient way to calculate the R-r in QCD or, turning the problem around, to find 
as{m1) from experimentally known R^- is to transform the integral in (1) to the integral 
over the contour in the complex s-plane going couterclockwise around the circle \s\ = 



10 — 13 : 
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The polarization operator is given by the sum of perturbative and nonperturbative terms. 
If we restrict ourselves by OPE terms, then 
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Consider at first the perturbative part. For its calculation it is convenient to use Adler 
function which is perturbatively constructed as an expansion in couphng constant 
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which is known up to the 4-loop term in MS renormalization scheme : Kq = Ki = 1 and 
K2 = 1.64 = 6.37 |15| for 3 flavors. The renormalization group equation for a{Q^) 

reads: 

'^'^ -/3(a) = -J]/3.a"+2 = - T""^^ ^ (9) 



dlnQ"^ 



n>0 



(3{a) 



In MS scheme for 3 flavors (3o = 9/4, /3i = 4, /Js = 10.06, (3^ = 47.23 0. This allows us 
to get the perturbative contribution to the polarization operator explicitly at any order of 
perturbation theory: 
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Figure 1: Real and imaginary parts of aMs('^)/^ exact numerical solution of RG equation 
(§) on real axes for different number of loops. The initial condition is chosen = 0.355 at 
s = — m^, Nf = 3. Vertical dotted lines display the position of the unphysical singularity at 
s = —Qq for each approximation (4-^1 from left to right). 



Let us put /i^ 



and choose some value a{m1). From (|^) we can find a(Q^) for any 



and by analytical continuation at any s. Computing the integral (|TUp it is possible 
to find the perturbative part of n(s) as a function of a{s) in the whole complex s-plane. 
The substitution of n(s) into (H) gives (up to the power corrections) the dependence of Rr 
It must be stressed, that in this calculation no expansion in inverse powers of 
is performed: only the validity of expansion series in (|) and (|) is assumed^. Such 
representation has a serious advantage: on the right semiaxes, i.e. in the physical region, 
there is no expansion in vr/ In (Q^/A^), which is not small at intermediate Q^. For instance 
in the next to leading order 



on a[m 
lnQ2 



27rlmn(s + i0) = 1 + 
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instead of 27rlmn(s + i0) = 1 



/?oln(s/A2) 

(Eq. ([TTl) was first obtained in 



which would follow in case of small ti / \n.{s/ h?). (Eq. (|11|) was first obtained in |]18[, the 
systematical method of analytical continuation from space-like to time-like region with sum- 
mation of TT^-terms was suggested in ||I9| and developed in pOf.) In the higher order, where 



a(s) cannot be expressed via ln(s/A2) in terms of elementary functions, this analysis is 
performed numerically. 

It is well known, that in the 1-loop approximation of /3 function the coupling a(Q^) 
has an infrared pole at some = Qq (in some conventions coinciding with A^). In the 

^Such way of calculation in |^ was called contour- improved fixed-order perturbation theory. 
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n-loop approximation (n > 1) instead of pole a branch cut appears with a singularity 
(1 — Q^/Qo)~^^"- "^^^ position of the singularity is given by 



da 



(12) 



Near the singularity the last term in the expansion of j3{a) (P) dominates and gives the 
aforementioned behavior. To illustrate the behavior of the running coupling constant, we 
plotted real and imaginary part of a^/Tr for n = 1,2,3,4-loop /3-function in Fig |I[ It 
demonstrates, that for real positive s the difference between various approximations is almost 
unnoticable beyond 2-nd loop and the expansion in inverse | ln(s/A^) — i7r\ works well. At 
the same time the behavior in the unphysical cut strongly depends on the number of loops 
and cannot be described by some simple approximation. Only at s < 1 GeV^ 2-4 loop 
calculations more or less coincide. 

Let us turn now to OPE terms in (|^. The contribution of the operators up to dimension 
8 have been computed theoretically: 
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The contribution of = 2 operator due to nonzero quark masses niu^ is negligible and 
omitted here. We have also neglected the D = A quark condensate 2(m„ + m^) {qq) which is 
an order of magnitude less than the gluonic condensate. The coefficients in front of D = 4, 6 

2T], hereafter cited as SVZ. The ag-correction to the 
-corrections to D = Q operator were calculated [p^ ; 



a. 



operators have been computed in 
D = A operator were found in p2 
ambiguities among them were also discussed there. 

Few comments about the operator Oq are in order. In nonfactorized form without 



corrections it looks as follows 211: 



2TiaA {u-i^X''d){d-f^X''u) + (M757A.A''d)(rf757MA"M) 



+ -{u^.X'^u + d^.X'-d) Y.iqi.X'^q) 



',,d,s 



(14) 



After factorization three terms in ([T^) give the following contributions: 



(Oe) = 4™^ {qq)' 



c 



(15) 



where Nc is the number of colors. SVZ assumed that the accuracy of the factorization 
procedure is of order A*""^ ~ 10% in case of V^-correlator, where the coefficient in the second 
brackets in ([T5|) is equal to —7/9. Remind that inV — A correlator the first term has opposite 
sign and the third term is absent, so the accuracy of the f actor ized operator 0^~^ is at least 
not worse, than in V case. On the other hand in the V + A correlator two comparatively 
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large terms cancel each other under the factorization assumption in (pISl). Consequently 
the accuracy of the formula (|T3p for the operator Oq is less, perhaps 20 — 30%. Large 
a^-corrections to all independent D = Q operators |2^ can only increase the errors. 

The numerical value oi D = Q operator can be estimated, for instance, with help of our 
previous analysis oiV — A sum rules : 

64 

(O^-^) = - — Tra, {qqf x 1.3 = - (6.8 ± 2.1) x 10"^ GeV^ (16) 
9 

The coefficient 1.3 stands for the a^-corrections. We find: 

(Oe) = (1.3 ±0.5) X 10-3 GeV^ , (17) 

The dimension 8 operators come from many different diagrams, which can be labeled by 
the number of quarks in vacuum. The purely gluonic condensates are suppressed by the loop 
factor ~ fts/TT and are neglected on this ground. The 4-quark operators, computed in ^ 
and 1^, vanish in the sum V + A after factorization. The uncertainty of this cancellation can 
be estimated as ~ 10% of O^""^, which is about 10^'^ GeV®. The 2-quark operators have the 
same sign in V and A correlators. They have been computed in (we have performed the 
calculation independently to confirm this result) and can be written in the following form: 



U 



+ ^ u{D^J)u - lu-f" [Ga^fs, J'^]u^+ {u^ d) (18) 

where = DyGi^^y = nagX'^ J^qi^'lf^'^'^l)- The last two terms can be factorized and brought 
to the form Tias{qq){qGq). However the leading in the number of colors A^^ terms cancel each 
other and only the terms ~ N'"^ are left. It has been shown in 0, that the factorization 
of D = 8 operators is not unambiguous at this level of accuracy. Taking the value of 
the operator {qGq) from |]2^, |[, we may estimate the upper limit of the operator (pTsP as 
KOg^'*)! < 10~'*GeV^, which is tiny. So, for the upper limit of the total D = 8 operator we 
shall use the estimation \ {Ps) \ < lO^^GeV^. 

It is worth mentioned, that the D = 6, 8 operators in + A polarization function are 
much smaller, than in or A separately. 

We are now in position to calculate as{Tn^) from the experiment. We take the most 
recent data on the total hadronic decay ratio Rr and the ratio of decays with odd number 
of strange mesons r~ X{S = —1)i't [H, HI]: 



= 3.636 ± 0.021 , i?^,^ = 0.161 ± 0.007 (19) 

In our analysis we subtract Rr^s to avoid the interference with additional parameters, in 
particular the mass of s-quark. One obtains 

Rry+A = 5\V^d\^SEw ( 1 + + + ^v+A ) + = 3.475 ± 0.022 (20) 

where AR^^ is given by (Kf). We use conventional in r-literature notations of fractional 



corrections 6. The electromagnetic correction is 6'^^^ = ■^aemi'mr) = 0.001 |Q, the D = Q 
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Figure 2: Perturbative fractional correction 6^^^ versus as{m'^) and as{Tn?z) in conventional 
and analytical approach in 3-loop approximation. In the width of the experimental strip the 
theoretical uncertainty of the operator (Og) is included. 

operator correction Sy\j^ = — (5 ± 2) x 10^^ as follows from our analysis, in agreement with 
the estimation obtained in [1^. From ( j20|) we separate out the perturbative correction: 

1 + ^(0) = 1.206 ±0.010 (21) 

All errors in here are added in quadratures (perhaps, such procedure underestimates the 
total error, may be by a factor 2). 

The calculation of asim^) corresponding to 5^^^ were performed according to the method 
described above. The dependence of 1 + 5^°^ on as{ml) (and on q;(M|), to compare with 
other data) for 3-loop /3-function and 3-loop Adler function is shown of Fig ^. It follows 
from Fig 0: 

asiml) = 0.355 ±0.025 (22) 

The estimation of the error in (^2]) was done with care. Because of asymptotic character 
of perturbative series (||) and (||) the higher loop contribution could be as large as the 
contribution of the last terms, namely K^a"^ and /32a^. They result to uncertainty 0.015 — 0.20 
in as{ml), depending on its central value. Taking into account the uncertainty (^) in 5^^\ 
we obtain the error in (|22|) . Furthermore, we have performed 2 and 4-loop calculation of 
as{m'^). The unknown 4-loop coefficient in Adler function (^ was taken equal to = 50 
(cf. its estimations |^). For each given the 4-loop as(m^) is by 0.005 lower than 3-loop 
value, while 2-loop ^^(m^) is by 0.02 higher. These results are within the error range ([2^). If 
some nonperturbative terms beyond OPE exist (e.g. instantons), they would also contribute 
to the error in (|22|) . In section 4 it will be shown, that the value ^^(m^) close to the lower 
limit of ( p2D satisfies sum rules at low much better. 
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as(m^) and analytical QCD 



Shirkov and Solovtsov forwarded the idea of analytical QCD. According to it the coupling 
constant as{Q'^) is calculated by renormalization group in the space-like region 

g2 > 0. Then, 

by analytical continuation to s = —Q^ > 0, as{s) was found, in particular its imaginary part 
Imas(s) on the right semiaxes. It was assumed, that as{s) is an analytical function in the 
complex s-plane with a cut along the right semiaxes < s < oo. The analytical c^s (s)an is 
then defined in the whole complex s-plane by dispersion relation: 

tt^Sjan = - / —ds (23) 

TT Jo S' 

Since the lower limit in this integral is put to zero, as(s)an indeed has no unphysical sin- 
gularities (poles, cuts etc) at > 0. The idea of analytical QCD has been developed in 



many papers, see e.g. |Q and for review 0]. In particular the calculations of a,,(m 



'an 



from r-decay data were performed within the framework of analytical QCD in 

A related approach was suggested by two of us (B.C. and B.I.) in We started 



from well-known theorem, that the polarization operator for e e~ annihilation n(s) is an 
analytical function of s in complex s-plane with a cut along positive semiaxes and assumed 
that these analytical properties take place separately for perturbative and nonperturbative 
parts of n(s). In the first order of this hypothesis is equivalent to analytical QCD while 
in higher orders it may be more general. 

Let us calculate as(^^)an in the framework of analytical QCD from the same experi- 
mental data, i.e. 6^^^ given by (0). The only (but important) difference from the previous 
calculation is the following. The coupling as(s)an is an analytical function of s with a cut 
from s = to s = oo. Consequently the contour integral in (H) is now equal to the orig- 
inal integral (^ with Imn(s) over real positive axes. In the previous calculation, if such 
transformation is performed, the integral would run from s = —Qq to m^. Qualitatively it 
leads to much smaller in analytical QCD than in conventional approach with the same 
«<i(m^), or vice versa, the same Rr corresponds to much larger as(s)an- Direct numerical 
calculation confirms this expectation. The dependence of 1 -|- S^°^ versus as(s)an is also 
displayed in Fig 0. It is seen, that in order to get experimental value of S^^^ in analytical 
QCD one should take as(m^)an ~ 1-5 — 2.0, which corresponds to as(m|) ^ 0.15 in strong 
contradiction with the world average as{m%) = 0.119 ± 0.002 (The previous calculation 
of r-decay performed with less certainty, demonstrated the same trend: in particular 
A^^) = 700 — 900 MeV, much larger, than in standard calculations.) 

It the recent paper an attempt was made to save the analytical QCD in case of 
vector polarization operator and to obtain the agreement with ALEPH data on vector Adler 
D-function by assuming large quark masses m„ = ma = 250 MeV and some form of Coulomb- 
like quark-antiquark interaction. This hypothesis, however, is in strong contradiction with 
all results following from well established partial conservation of axial current (PCAC) and 
chiral theory. For example, Gell-Mann-Oakes-Renner relation and K/n mass ratios would 
be violated by an order of magnitude, Goldberg- Treitman relation cannot be proved etc. 
Also many sum rules for V — A polarization operator would disagree with the data. 



Therefore we come to the conclusion, that analytical QCD in any form, or ||36| is in 
strong contradiction with experiment and must be abandoned. 
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Figure 3: Region of validity of the perturbation theory (PT) and operator product expansion 
(OPE) 

4 Check of QCD at low for V + A correlators by 
using the sum rules 

Let us turn now to study of the ^ + A correlator in the domain of low Q^, where the OPE 
terms play much more essential role, than in the determination of i?,-. A general remark is in 
order here. As was mentioned in and stressed recently by Shifman the condensates 
cannot be defined in rigorous way, because there is some arbitrariness in the separation of 
their contributions from perturbative part. Usually |]3S|, ^ they are defined by introduction 
of some normalization point /i^ with the magnitude of few A^. The integration over momenta 
in the domain below /i^ is addressed to condensates, above //^ — to perturbation theory. 
In such formulation the condensates are /i-dependent {Od) = {Od) ^ and, strictly speaking, 
they also depend on the way how the infrared cut-off /i^ is introduced. The problem becomes 
more severe when the perturbative expansion is performed up to higher order terms and the 
calculation pretends on high precision. Mention, that this remark does not refer to chirality 
violating condensates, because perturbative terms do not contribute to chirality violating 
structures. For this reason, in principle, chirality violating condensates, e.g. (0|gg|0), can be 
determined with higher precission, than chirality conserving ones. Here we use the definition 
of condensates, which can be called n-loop condensates. As was formulated in Section 2, we 
treat the renormalization group equation (^ and the equation for polarization operator ( pUj ) 
in n-loop approximation as exact ones; the expansion in inverse logarithms is not performed. 
Specific values of condensates are referred to such procedure. Of course, their numerical 
values depend on the accounted number of loops; that is why the condensates, defined in 
this way, are called n-loop condensates. 

Consider the polarization operator 11 = H^^^ + H^'', defined in (|I|) and its imaginary 
part 

uj{s) = vi{s) + ai{s) + ao{s) = 27rlmn(s + i0) (24) 
1 at s ^ oo. Any sum rule can be written in the following form: 



In parton model w{s) 



/so r 
fis)ujc^p{s)ds = in j) f (s) Utheovis) ds (25) 

where f{s) is some analytical in the integration region function. In what follows we use 
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uJcxp{s), obtained from r-decay invariant mass spectra published in [||] for < s < with 
step ds = 0.05 GeV^. The experimental error of the integral (U^ is computed as the double 



integral with the covariance matrix uj{s)uj{s') — iJ{s)uJ{s'), which also can be obtained from 
the data available in [|l| . In the theoretical integral in (^) the contour goes from sq + iO to 
sq — iO counterclockwise around all poles and cuts of theoretical correlator n(s), see Fig ^ 
Because of Cauchy theorem the unphysical cut must be inside the integration contour. 

The choice of the function f{s) in (pSj) is actually a matter of taste. At first let us 
consider usual Borel transformation: 

= f' e-'"-".,^,{s) ^ = B„.(M^) + 2.' ^2 (26) 

JO n ^ ' 

We separated out the purely perturbative contribution i?pt, which is computed numerically 
according to ( pS]) and (p|-p!0|). Remind that Borel transformation improves the convergence 
of OPE series because of the factors — 1)! in front of operators and suppresses the 
contribution of high-energy tail, where the experimental error is large. But it does not 
suppress the unphysical perturbative cut, the main source of the error in this approach, even 
increase it since e~^l^ > 1 for s < 0. So the perturbative part i?pt(M^) can be reliably 
calculated only for 0.8 — 1 GeV^ and higher; below this value the influence of the 

unphysical cut is out of control. 

Both 5exp and Spt in 3-loop approximation for a<j(m^) = 0.355 and 0.330 are shown 
m Fig |. The shaded areas display the theoretical error. They are taken equal to the 
contribution of the last term in the perturbative Adler function expansion K^a^ We 
have also performed the calculation with 4-loop /5-function and ^^^4 = 50 ± 50, but the result 
is very close to the 3-loop one, since positive contribution of the term K^a^ compensates 
small decrease in the coupling a. Since this result is observed by us in many other sum rules, 
we shall not give the 4-loop calculations later on and estimate the theoretical error for any 
given a(m^) as the contribution of K^a^. 

As follows from the analysis in Section 2, for > 1 GeV^ the contribution of D = 6, 8 
operators to the Borel transform (^) is small in K + A channel, while the contribution of 
the D = 4 condensate must be positive (we assume Os-corrections included in the operators 



(02n) in ( 126|) and later). So the theoretical curve must go below experimental one. The 
result shown in Fig |^ is in favor of lower value of the coupling constant a;s(m^) = 0.33. 
Literally the theoretical curve (perturbative at as(m^) = 0.33 plus the contribution of 
and Oq operators) agrees with experiment starting from = 1.1 GeV^. If the uncertainties 
in perturbative contributions are taken into account (shaded area in Fig ^ the agreement 
may start earlier, at = 1 GeV^. 

The Borel transformation on Fig |^ includes the contributions of different operators. Al- 
though it is difficult to separate the perturbative part from the OPE one, the contributions 
of different operators can be separated from each other. One way is to differentiate the Borel 
transformation by M^. This however leads to the certain loss in the accuracy of the experi- 
mental integral, since the growing power term ~ s" appears in the integral. So we apply the 
method used in Q for V — A sum rules, namely the Borel transformation in complex 
plane. 

Let us consider the Borel transform i?(M^) (^) at some complex = M^e'"^, < 
< 7r/2. If the phase (j) is taken close to 7r/2 then the contribution of the high-energy 
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Figure 4: Borel transformation (|6|) B^^p{M^) and Bpt{M^) for ^.(m^) = 0.355 and 0.330. 
The dashed hne displays the OPE contribution added to the 0.330-perturbative curve. The 
contribution of the operators D = 4 (standard SVZ value) and D = 6 (central value of (p!7|)) 
with respect to 1 are shown separately in the box. 



tail becomes high. So we restrict ourselves by the values < 7r/4 for the exponent to be 
decreasing enough. The real part of the Borel transform at = it/ 6 does not contain the 
D = 6 operator: 

ReB,^p{M'e'-/') = ReBp,{M'e'-/') + vr^ ^ (27) 

The contribution of (Og) is less than 0.5% to the perturbative term and neglected here. The 
results are shown in Fig |^. Again it is still difficult to accommodate positive value of the 
gluonic condensate to the coupling ^^(m^) = 0.355 and higher. If we accept the lower value 
of asim"^), we get the following restriction on the value of the gluonic condensate: 

^G'l^vG'Jl^) = 0.006 ±0.012 GeV\ ^.(m^) = 0.330 and M'^>0.8GeV^ (28) 



The theoretical and experimental errors are added together in (28 



The real part of the Borel transform at = 7r/4 does not contain the D = A operator: 

ReB,^p{M'e^-/') = ReBp,{M\^^'') - (29) 

The results are shown in Fig §b. The perturbative curve at as = 0.330 is below the data. 
If we would take this curve as an exact one, without accounting the perturbative errors, 
then from (^) we would conclude, that (Oq) < 0, which in some contradiction with (|13|, p!7| ). 
However the account for the perturbative errors makes the situation different, but uncertain. 
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Figure 5: Real part of the Borel transform (pGf) along the rays at the angles = 7r/6 and 7r/4 
to the real axes. The dash line corresponds to the gluonic condensate given by the central 
value of (H). 



Since the value of the (Oe) contribution to ( p9D is very small 

^2 (Oe^ = (0.9 ± 0.4) X 10-2 ^ (30) 

then by accounting the perturbative errors it is possible to satisfy the sum rule (|^) at 
positive {Oq) starting from > 0.8 GeV^. (In the narrow region near = 0.9 GeV^ the 
theoretical curve goes out of the data on 1.5 — 2 experimental error, but we do not consider 
this as a serious contradiction.) Unfortunately no definite conclusion about the value of Oq 
can be done from the Fig Qo. The only statement is that its value cannot exceed ([T7|) and 
probably is on the lower border of error. 



5 Correlator of vector currents 

Previously we considered the V+A correlators where the power corrections are small. Instead 
one could take pure vector current, (vector spectral function was published by ALEPH in 
40[|). This doesn't give us any new information with the r-decay data, since V—A correlators 



have already been analyzed in [Q. Moreover the accuracy of the vector current spectral 
function is less, than + A, since both currents are mixed in some channels with ii'-mesons 
and the number of events is twice less. 

However the analysis of the vector current correlator is important since it can also be 
performed with the experimental data on e~^e~ annihilation. The imaginary part of the 
electromagnetic current correlator, measured here, is related to the charged current correlator 
(0) by the isotopic symmetry. The statistical error in e~^e~ experiments is less than in r- 
decays because of significantly larger number of events. So it would be interesting to perform 
similar analysis with e~^e~ data, which is a matter for separate research. 
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Figure 6: Borel transformation for vector currents. 
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At first we consider usual Borel transformation for vector current correlator, since it was 
originally applied in ^ for the sum rule analysis. It is defined as ( PBD with the experimental 
spectral function cUgxp = 2t>i instead of Vi + ai + gq (the normalization is fi(s) — ^ 1/2 at 
s — > oo in parton model). Respectively, in the r.h.s. one should take the vector operators 
20^ = 0^+^ + O^-^, all O^-^ with D < 8 can be found in §. Numerical results are 
shown in Fig ^ The perturbative theoretical curves are the same as in Fig | with V + A 
correlator. The dashed lines display the contributions of the gluonic condensate given by 



20 



-5.5 X 10"' GeV 



and 20l 



qV- 



7 X 10"'^ GeV* added to the 0.330- 



perturbative curve. The contribution of each condensate is shown in the box below. Notice, 
that for such condensate values the total OPE contribution is small, since positive O4 and 
Og compensate negative Og. The agreement is observed for > 0.8 GeV^. 

Now we apply the method of Borel transformation along the rays to the vector polariza- 
tion operator to separate the contribution of different operators from each other. The Z) = 8 
operator is important here, so we shall separate 04^6 from Og. 

Borel transformation at low exponentially suppresses the contribution of large s 
domain, where the experimental error is high. Besides of this we may use the oscillating 
behavior of the complex exponent to further suppress the high-error points near s = m^. This 
would allow us to go to higher M^. Here the real part of B{M'^) has obvious advantage since 
the function cos {4> + cos0) has zero at s = and ~ 7r/4 already at ^ 1 GeV^, 

r2\ /'A _|_ _|_ COS0) in imaginary part is twice lower. So let 

0,7r/6,7r/4. Solving the system of linear equations. 



while the largest (in M^) zero of sin 
us take three different angles, say, 
we get: 



Re 



5(0) 



.2- 
Re 



-5(0) + 2B{n/6) 



B{tt/A) 

v^- 1 
2 5(7r/4) 



V2- 1 



p.t. 
p.t. 



2r 



9±. 

2M6 



(31) 
(32) 



For brevity we write here -8(0) instead of i?cxp(^^e*'^), "p.t." stands for perturbative con- 
tribution. The results for the equations (PT| , P^]) are shown in Figs |^,b respectively. 

Fig 1^ demonstrates, that the vector sum rule is satisfied at ^^(m^) = 0.330 and gluonic 
condensate ( pSj) (although higher values of gluonic condensate, e.g. SVZ value still do not 
contradict the data). Fig ^b shows, that the contribution works in right direction: 
its addition to 0.330-perturbative curve shrinks the disagreement between the theory and 
experiment. However, some discrepancy (about 0.04, i.e. 0.1% in the worst case) still persist. 
It may be addressed either to the uncertainty in ^^(m^) — a slightly higher value would be 
desirable, or to the underestimation of (in absolute value: is negative) by 20 — 30%, or 
both. Remind, that the numerical values of condensates depend on the way, how the infrared 
region is treated (Og is chirality conserving). We are considering here 3- loop condensates, 
defined in Sec. 4. The 0\ value was taken equal to 7/18 of obtained from V — A 

data analysis [Q, where perturbative terms are absent, and some difference is not excluded. 
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Figure 8: Sum rule ( pSj) with O2 but without 04,6 



6 The check for renormalon-type terms 

In the asymptotic perturbative series a special part of terms — renormalons (infrared and 
ultraviolet) is often separated and the summation of them is performed (for a recent review 
see In such sum the term appears proportional to l/Q^ at large and looking as a 

contribution of -D = 2 operator. (In OPE the D = 2 operator is proportional to and is 
very small.) Renormalons conserve chirality and may contribute ioV + A but not ioV — A. 
Unfortunately, the coefficient in front of the term of the renormalon origin cannot be 



calculated reliably. (In |^ it was claimed, that the renormalons are totally absent in the 
perturbative series asymptotics and therefore this coefficient is zero.) In recent paper [43| 
the hypothesis was suggested, that infrared renormalons result in substitution 



- ^ - 1 - 1-0575^ 



(33) 



in the first correction to polarization operator or Adler function (the Q^-dependence of 
as was not accounted in [Q). In (|33| ) is tachionic gluon mass, A^ < and for its value 
the estimation was found: 

- A^ = (0.2 - 0.5) GeV^ (34) 
The authors of could not discriminate even the highest value A^ = — O.SGeV^. 



Let us try to find the restriction on O2 operator from the sum rule for V + A correlator 
in the complex g^-plane from ALEPH data. (We call it for brevity O2, although it is not the 
D = 2 operator which stands in OPE.) As we did in previous section, for this purpose we 
take the real part of the Borel transform (|26| ) i?(M^e*'^) at the angles = 0,7r/6,7r/4 and 
separate the operator O2 from O^q: 



Re 



B{0) - 2B{-n/Q) + ^/2B{-k/A) 
2- V3 



p.t. 



2ti' 



M2 



(35) 
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The experimental and perturbative parts of this combination are plotted in Fig p. 

The sum rule (H^) shown in Fig R gives the following value of the dimension 2 operator: 



O2 = (1.0 ± 1.5) X IQ-^ GeV^ , ^.(m^) = 0.33 (36) 

We got this estimation at = 1 GeV^, where experimental error is minimal. In the model 
of g 



a 



O2 = -1.05 — A' (37) 



TT 



At as(l GeV )/vr = 0.18, corresponding to ^^(m^) = 0.33, there follows the restriction from 

- A^ = (5 ± 8) X 10-2 GeV^ (38) 



which is few times smaller than even the lower limit in (^^. Notice, that similar restrictions 
on the value oi D = 2 operator have been obtained in E3| from other sum rules. 



7 Instanton corrections 



Some nonperturbative features of QCD may be described in so called instanton gas model 
(see |^5| for extensive review and the collection of related papers in [0). Namely, one 
computes the correlators in the S'f/(2)-instanton field embedded in the SU{3) color group. 
In particular, the 2-point correlator of the vector currents has been computed long ago [0. 
Apart from usual tree-level correlator ~ InQ^ has a correction which depends on the 
instanton position and radius p. In the instanton gas model these parameters are integrated 
out. The radius is averaged over some concentration n{p), for which one or another model 
is used. Concerning the 2-point correlator of charged axial currents, the only difference 
from the vector case is that the term with 0-modes must be taken with opposite sign. In 
coordinate representation the answer can be expressed in terms of elementary functions, see 
An attempt to compare the instanton correlators with ALEPH data in coordinate space 
has been undertaken in |4^ . 

We shall work in momentum space. Here the instanton correction to the spin-J parts 
n^"^'' of the correlator (|l|) can be written in the following form: 



n?,Lt(?') 
n2U(g^) 

(1) ^^2^ 



n 



A, inst 



dpn{p) 



dpn{p) 



+ V^p'Gll -pV 



3g4 

4 4p2 



— — 2^Apy~^ 



(0) 



■'■■'■y.mstW 



1/2 
0,0,-2 







(39) 



Here is modified Bessel function, G^Jyz\ 



is Meijer function. Definitions, properties 

In particular the 

Of) can be written as the following series: 



and approximations of Meijer functions can be found, for instance, in ^9 
function in 



0FG?° I z 



1/2 
0,0,-2 



4 2 

32;2 z 



2^ 



fc=0 



r(A; + l/2) 



r2(A; + 1) r(A; + 3) 



X 
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x| [In^ + '4){k + 1/2) - 2^(fc + 1) - ^{k + ?,y 
+ ij'{k + 1/2) - 2ij'{k + 1) - ij'{k + 3) 



(40) 



^30 



\z\ > 1 



(41) 



where ip{z) = r'{z)/r{z). For large \z\ one can obtain its approximation by the saddle-point 
method: 

1/2 ■ 
0,0,-2 

The formulas (|39D should be treated in the following way. One adds Ilinst to usual polariza- 
tion operator (^ with perturbative and OPE terms. But the terms ~ 1/g^ must be absorbed 
by the operator O4 in (|^), since the gluonic condensate {G^) is averaged over all field config- 
urations, including the instanton one. Notice negative sign before 1/g^ in (|39|) . It happens 
because the negative contribution of the quark condensate {mqq) in the instanton field ex- 
ceeds positive contribution of the gluonic condensate (G^). In real world {mqq) is negligible 
at g2 ^ 1 GeV^. 

The correlators (|39|) possess appropriate analytical properties, they have a cut along 
positive real axes: 



Snst(9'+^0) 



imn^^^ '-^ 



dpn{p)n'/'p'GlUp'q 



1/2 
0,0,-2 



Imn?) (g^ + ^O) 



dpn{p) 



27rV' 



(42) 
(43) 



We shall consider below the instanton concentration advocated by Shuryak (see |^5| and 
references therein). It is a model with fixed instanton radius (RILM model in RH]): 



n(p) = noS{p- po) 



(44) 



From we take the numbers: 

Po = 1/3 fm = 1.7 GeV"^ 



no 



1 fm 



1.5 X 10"^GeV^ 



(45) 



Now we consider the instanton contribution to the r-decay branching ratio (§). Since 
the instanton correlator (|39| ) has 1/g^ singular term in the expansion near (see (|iO|)), the 
integrals must be taken over the circle, like in (^. In the instanton model the function ao{s) 
differs from experimental (5-function, which gives small correction (^. So we shall ignore 
the last term in (Q) and consider the integral with Hy^^ -|- II^'* in (^. Here we need the 
following formulas for the circle integrals, which can be rigorously obtained from the series 
representation of the Meijer function 



2^ 



-f-l Git 



\s\=so ^0 

i 



2tt 



'-^13 

|s|=so So So 
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2^ 



\s\=so So 



^S 

^13 





1/2 \ 
0,0,-2y' 




1/2 \ 


0,0,-2;* 




1/2 \ 


0,0,-2;* 



Gii [^p'so 
4 

-Gil (p'so 



-k, 1/2 
0,0,-2,-A;- 1 



k>2 



+ ^24 yP sq 

1/2 
0,-1,-2 



-1,1/2 
0,0,-2,-2 



(46) 
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Figure 9: The instanton correction to the r decay ratio versus po (a) and "versus r mass" 
(b) for no = 1.5 X lO'^ GeV^ 



The first term in the r.h.s. of the second equation looks hke the contribution of D = 4 
operator, but in fact it is not. Indeed, all expressions in the r.h.s. of ( ^61) have the same LO 
term of the asymptotic expansion for large sq, equal to — sin (2py^)/(y^p^SQ). However 
for k > 3 the accuracy of this approximation is bad and exact values of Meijer functions 
should be used for numerical evaluations. 

With help of (^) the instanton correction to the r-decay branching ratio can be brought 
to the following form: 



48 TT^/^ 



1/2 
0,-1,-4 



sin {2pQm^ 



(47) 



Since the parameters (^5]) are determined quite approximately, we may explore the depen- 
dence of ^inst on them. The 5inst versus po for fixed Uq (^) is shown in Fig 

As seen from Fig |^a the instanton correction to hadronic r-decay is extremely small 
except for unreliably low value of the instanton radius po < 1.5GeV~^. At the favorable 
value 1^ Po = 1.7GeV~^ the instanton correction to is almost exactly zero. (Of course, 
smaller values of no than (|45| ) are also allowed.) This fact confirms our calculations of a;s(m^) 
(Sec. 2), where the instanton corrections were not taken into account. 

The result (^^ can be used in another way. Namely, the r mass can be considered as 
free parameter sq- The dependence of the fractional corrections 5'-°'' and 5^330 + "^mst on 
So is shown in Fig ^30. The result strongly depends on the instanton radius and rather 
essentially on the density no- For po = 1.7GeV~^ and no = lfm~^ (^5]) the instanton curve 
is outside the errors already at Sq ~ 2GeV^, where the perturbation theory is expected to 



^The Fig ^ can be compared with Figure 15 in ALEPH paper The discrepancy between theoretical 
curves at sq < 1 GeV^ is explained by different approximations: we used 3- loop perturbation theory, while 
the authors of |Q] used 4-loop one with K4 — 50 ± 50. 
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work. Therefore Fig Pd shows, that in RILM model the instanton radius must be larger 
(say, po = 2.5 GeV~^) or the instanton density much (2 — 3 times) lower. The contribution 
of D = 6 operator Sy}^j^ is not shown on Fig ^b. It is equal Sy]_^ = — (5 ± 2) x 10~'^(m^/so)^ 
and quite large at sq < 1.5 GeV^. 

Consequently in this approach the perturbation theory + OPE + RILM (at not very large 
Po) cannot satisfactory describe the data at Sq < 1.5 GeV^. Since the instanton contribution 
is large here, we disbelieve all the results, obtained by the method of variable r-mass in this 
domain. (Perhaps, the shadowed region in Fig ^ is of importance in this method at low Sq.) 

The r decay ratio is not sensitive to the gluonic condensate. Let us consider now the 
sum rules which depend on it. The Borel transformation of the instanton part is: 



27Ti 



An' 



-s/M^jt{1) (^\_^ 



oo 



dpn{p) 



4 



3M4 



1/2 
0,-2 



(4J 



The integration contour goes around the cut from s = 
term ~ here comes from the term ~ 1/g^ in (^ 



l-oo + iO to s = +00 — iO. The 
it must be included in the (O4) 



contribution in (^61). The Meijer function in (^Sf) has the asymptotics 
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1/2 
0,-2 



\z\ > 1 



and strongly suppressed at > 0.8 GeV^. We calculated the instanton contribution to all 
Borel-like sum rules used here; it is indeed negligible compared to the errors. Consequently 
the results of previous sections remain unchanged. 

However the spectral moments sum rules, often used in r-decay data analysis [Q, can 
be quite sensitive to the instanton corrections. Let us consider the following sum rule, 
constructed in this way: 



4 r^A 

Jo Sq 



i-£: 

si 



00, 



exp I 



p.t. _ Stt^^ + 167r2 

^^0 



POO 



dpn{p) p"" 



3p^sl 



24 



24 



-3,1/2 



0,0,-2,-4 



(49) 



-1,1/2 
0,0,-2,-2 

The integral ( ^Dj ) is normalized to 1 in parton model. It does not depend on the D = 6 
operator, the factor 1 — s^/ Sq is introduced to suppress large experimental errors for large Sq. 
Remind our convention: the contribution of the term ~ 1/g'^ in Ilinst ( |5D| ) is included into 
the operator (O4) in (^9]). The contribution of different parts of eq. ( ^9D are plotted versus 
So in Fig |10|. Since the weight function in the integral vanishes at s = 0, the contribution 
of unphysical cut is suppressed. So the theoretical errors are diminished here as well as 
the sensitivity on various perturbative parameters. The theoretical curve is shown as single 
shaded area, which includes both the uncertainty of Q;s(m^) and the error ±^30^ for each 
asiml). 

The operator (O4) enters with negative sign in (^), so the theoretical curve must go 
above experimental one. This is certainly not the case if the instanton corrections are not 
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Figure 10: Sum rule (|i9|). Experimental, pure perturbative and "perturbative + instanton" 
parts are shown. The O4 contribution is not taken into account. 

taken into account. For po = 2.1 GeV^ the theoretical and experimental results are in 
good agreement for (O4) = 0. By increasing the instanton density no, positive values of 
(O4) become possible. In this aspect sum rules (^) with varying and (|49|) are not in 
agreement: (|47| ) favors small uq while (^9]) prefers large uq. 

These results are, however, not convincing. The main conclusion, coming from consider- 
ation of spectral moments sum rules, is that they are not suitable for QCD analysis untill 
we have a complete theory. (This statement surely refers also to the method, where r-mass 
is considered as free parameter.) The same situation took place for V — A correlators: the 
spectral moments sum rules worked only at the circle radius sq > 2 GeV^ . 

8 Conclusion 

The goal of this paper was to confront the recent precise experimental data on hadronic r- 
decay with QCD calculations at low and to check the basic aspects of QCD: perturbative 
series, OPE as well as various nonperturbative QCD approaches. The data present the 
imaginary part of polarization operators lmIlv^A{s) , s = at < s < m^. If some procedure 
is applied to suppress or nullify the influence of high energy domain (Borel transformation, 
integration over closed circle in complex s-plane), then with help of dispersion relation the 
values of IIv,a{s) is the whole complex s-plane at low \s\ can be found from experiment. (By 
low |s| we mean |s| < 2 — 3 GeV^.) These experimental values of Ilv,Ais) can be compared 
with theoretical calculations in the domain of complex s-plane, where QCD describes the 
data well enough, in order to find the values of QCD parameters: and condensates. 

In this program was realized for IIv-a polarization operator and the values of dimen- 
sion 6 and 8 condensates were found. In this paper Ily+A and IIv polarization operators were 
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studied, where perturbative contribution is dominant (unlike Ily.^ which is given entirely 
by condensates). It must be stressed, that the present situation has changed drastically in 
comparison with earlier study of similar problem In |^ the perturbative contribution was 
much less essential and the authors could restrict theirselves to LO term only. In this paper 
the perturbative calculations were performed in 3 and 4 loop approximation. The unphysical 
cut in the complex s-plane in perturbative part of the polarization operator was taken into 
account and the calculations (at least partly) were performed in such a way, which allows 
one to minimize its influence (e.g. the Borel transformation along the rays, going from the 
origin at some angle). The terms of OPE were accounted up to dimension D = 8. It was 
shown that D = 8 contribution is very small in case of V + A correlator. The coincidence 
of theoretical and experimental values with accuracy better than 2% was required. Let us 
remind that usually the accuracy of standard QCD sum rule calculations is of order 10 — 15%. 
The following results have been obtained: 

1. The value of QCD coupling constant as{'m'^) = 0.355 ±0.025 was found from hadronic 
branching ratio Rr- It was shown, the sum rules at low |s| favor the value close to the 
lower error edge as(m^) = 0.330 corresponding to as{m\) = 0.118. 

2. It was demonstrated that QCD with inclusion of OPE terms is in agreement with 
the data at the values of complex Borel parameter |M^| > 0.8 — 1.0 GeV^ in the left 
complex half-plane. 

3. The restriction on the value of the gluonic condensate was found (^G^) = (0.006 ± 
0.012) GeV^ in comparison with standard SVZ value 0.012 GeV^. 

4. The value oi D = Q condensate found in P] is in agreement with V + A and V sum 
rules, but cannot be specified. 



5. The analytical perturbative QCD [|32|, ^ was compared with the data and it was 
demonstrated that this approach is in strong contradiction with experimental value of 
R-f- 

6. The restrictions onl/Q^ term in polarization operator of renormalon origin were found. 



much stronger, than in the previous investigation |43 



The instanton contributions to polarization operator were analyzed and compared with 
the data in the framework of the random instanton liquid model (RILM) It was 



shown that the instanton contribution to R-r is very small, the same is true for Borel 
sum rules. However their contributions can be significant to the spectral moments sum 
rules, often used in r-decay data analysis. 

8. It was found that the method of spectral moments (integration over the circle with 
a polynomial) is less effective in the study of the polarization operators at low Q^, 
than Borel sum rule because of larger contribution not given by OPE nonperturbative 
corrections (see Sec. 7 and 0). 

We believe, that the results of this paper will serve for improving the QCD sum rules 
method. 
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